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Abstract 

We study a system composed of fermions (electrons), hopping on a square lattice, 
and of immobile particles (ions), that is described by the spinless Falicov-Kimball 
Hamiltonian augmented by a next-nearest-neighbor attractive interaction between 
the ions (a nearest-neighbor repulsive interaction between the ions can be included 
and does not alter the results). A part of the grand-canonical phase diagram of this 
system is constructed rigorously, when the coupling between the electrons and ions 
is much stronger than the hopping intensity of electrons. The obtained diagram 
implies that, at least for a few rational densities of particles, by increasing the 
hopping intensity the system can be driven from a state of phase separation to a 
state with a long-range order. This kind of transitions occurs also, when the hopping 
fermions are replaced by hopping hard-core bosons. 
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1 Introduction 

Since forty years, it is the Hubbard or the extended Hubbard models that 
are most frequently studied when properties of strongly correlated electrons 
are to be investigated. The central issue is the phase diagram of these mod- 
els. Despite the apparent simplicity of both models and concerted efforts of 
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many researchers, the complicated structure of these phase diagrams, even 
at zero temperature, has not been revealed completely and unquestionably; 
rigorous results are scarce. A few years ago Nakamura [1] predicted a new 
kind of quantum phase transition in the ground state of a half-filled extended 
Hubbard model chain. On a line, in the space of the two relevant interaction 
parameters (representing on-site and nearest-neighbor (n.n.) repulsion due to 
Coulomb forces, expressed in the units of the hopping intensity), for large 
values of these parameters, the system is in the phase-separated state, which 
is a mixture of a charge-density-wave phase and a spin-density-wave phase. 
The transition occurs when the parameters are decreased along the line, at a 
critical point where the parameters assume intermediate values (of the order 
of one). Thus, the transition region is hardly accessible by perturbation meth- 
ods. A new state exhibits a long-range order, the so called bond-order- wave. 
There are many aspects of the phase diagram to be studied in connection 
with the new conjectured transition, and they are vigorously discussed in the 
physics literature recently [2,3,4]. In this paper we address just one aspect of 
the new transition that we find remarkable. Namely, by increasing the hopping 
intensity of electrons the system is driven from a phase-separated state to a 
crystalline state with a long-range order. Our aim is to construct a model of 
interacting electrons (by the way we consider an analogous model of inter- 
acting hard-core bosons), where the existence of a transition exhibiting this 
feature can be demonstrated rigorously. The model to be studied is a simpli- 
fied version of the one band, spin 1/2 Hubbard model, known as the static 
approximation (one sort of electrons hops while the other sort is immobile), 
with the Hamiltonian H FK , extended by a next-nearest-neighbor (n.n.n.) at- 
tractive interaction between the immobile particles, given by the Hamiltonian 
V. Thus the total Hamiltonian of the system reads: 



H 



FK 



Ho = H FK + V, 
~t 2 ( c t c y + h - c ) + U Y, i c t c x - 1) s x , 



V = -— J2 s xS y , t>0, e>0. 
{x,y)2 



(1) 
(2) 

(3) 



In the above formulae, the underlying lattice is a square lattice, denoted A, 
consisting of sites x,y, . . . whose number is |A|, having the shape of a \J\A\ x 

\J\A\ torus. In (2,3) and below, the sums J2( x , y )^ i — 1)2,3, stand for the 
summation over all the i-th n.n. pairs of lattice sites, with each pair counted 
once. 

The subsystem of mobile spinless electrons is described in terms of creation 
and annihilation operators of an electron at site x: c+, c x , respectively, satis- 
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fying the canonical anticommutation relations, with t being the n.n. hopping 
intensity. The total electron-number operator is iV e = J2 X c t c x, and (with a 
little abuse of notation) the corresponding electron density is p e = iV e /|A|. 

The subsystem of ions is described by a collection of pseudo-spins {s x } xeA , 
with s x = 1, —1 (1 if the site x is occupied by an ion and —1 if it is empty), 
called the ion configurations. The total number of ions is Ni = J2xi s x + l)/2 
and the ion density is pi = A^/|A|. In contradistinction to the electron subsys- 
tem, the ions interact directly: two ions that occupy two n.n.n. sites attract 
each other, contributing the energy — e/16, with e > 0. 

Clearly, in the composite system, whose Hamiltonian is given by (1), with 
arbitrary electron- ion coupling U, the particle- number operators N e , N iy and 
pseudo-spins s x , are conserved. Therefore the description of the classical sub- 
system in terms of the ion configurations S = {s x } xeA remains valid. Whenever 
periodic configurations of pseudo-spins are considered it is assumed that A is 
sufficiently large, so that it accommodates an integer number of elementary 
cells. 

Nowadays, Hfk is widely known as the Hamiltonian of the spinless Falicov- 
Kimball model, a simplified version of the Hamiltonian put forward in [5]. A 
lot of results, including rigorous ones, like a proof of the existence of a phase 
transition [6,7] for instance, have been obtained for the system described by 
this Hamiltonian (a review and an extensive list of references can be found in 
[8,9]). 

In what follows, we shall study the ground-state phase diagram of the system 
defined by (1) in the grand-canonical ensemble. That is, let 

H (p e , p^ = H - p e N e - piNi, (4) 

where p e , pi are the chemical potentials of the electrons and ions, respectively, 
and let Es (p e , fa) be the ground-state energy of H (p e , p^), for a given config- 
uration S of the ions. Then, the ground-state energy of H (p e , pi), Eg (p e , Pi), 
is defined as E G (p e , p^) = mm{E s (p e , Pi) '■ S}. The minimum is attained at 
the set G of the ground-state configurations of ions. We shall determine the 
subsets of the (p e , /Xj)-plane, where G consists of periodic configurations of 
ions, uniformly in the size of the underlying square lattice. 

The paper is organized as follows. In the first subsection of the next sec- 
tion we provide a number of symmetry properties of the ground-state en- 
ergy, Es (p e , Pi), which facilitate the studies of the phase diagram that follow, 
and give the strong-coupling expansion of Es(p e ,Pi) in the case of hopping 
fermions (subsection 2.1). Then, in subsection 2.2 we construct the grand- 
canonical phase diagram due to the ground-state energy truncated at fourth 
order of the strong-coupling expansion in the case of hopping fermions (this is 
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sufficient to draw our conclusions). After that, in subsection 2.3 we consider the 
grand-canonical phase diagram due to the ground-state energy truncated at 
fourth order of the strong-coupling expansion in the case of hopping hard-core 
bosons. Finally, in Section 3 we discuss the implications of the phase diagram 
obtained with the truncated effective interaction, when the remainder of the 
strong-coupling expansion is taken into account. 



2 Grand-canonical phase diagram in the strong coupling regime 

2.1 Properties of E s (fa, fa) and its strong- coupling expansion 

In studies of grand-canonical phase diagrams an important role is played by 
unitary transformations that exchange particles and holes: c^c x — > 1 — c+c x 
and s x — > — s x , and for some (/i®, fa ) leave the Hamiltonian H (fa,, fa) invari- 
ant. For the electron subsystem such a role is played by the transformations: 
c+ — > e x c x , with e x — 1 at the even sublattice of A and e x — — 1 at the odd 
one. Clearly, since Hq is invariant under the joint hole-particle transforma- 
tion of electrons and ions, H (fa,, fa) is hole-particle invariant at the point 
(0,0). At the hole-particle symmetry point, the system under consideration 
has very special properties, which simplify studies of its phase diagram [7]. 
Moreover, by means of the defined above hole-particle transformations one 
can determine a number of symmetries of the grand-canonical phase diagram 
[10]. The peculiarity of the model is that the case of attraction (U < 0) and 
the case of repulsion (U > 0) are related by a unitary transformation (the 
hole-particle transformation for ions): if S is a ground-state configuration at 
(fa, fa) for U > 0, then —S is the ground-state configuration at (fa, —fa) f° r 
U < 0. Consequently, without any loss of generality one can fix the sign of 
the coupling constant U. Moreover (with the sign of U fixed), there is the 
inversion symmetry of the grand-canonical phase diagram, that is, if S is a 
ground-state configuration at (fie, fa), then — S is the ground-state configura- 
tion at (—fa, —fa)- Therefore, it is enough to determine the phase diagram in 
the half-plane specified by fixing the sign of one of the chemical potentials. 
Additional properties emerge in the strong-coupling regime, i.e. for \U\ > At 
[10]. Specifically, if the electron chemical potential is in the open interval 
\fa\ < \U\ — At, then for any two ion configurations S, S' the energy difference 
Es (fa, fa) — Es> (fa, fa) is constant along the lines fa±fa = const, with minus 
sign referring to the case of positive U (repulsion) while the plus sign — to 
the case of negative U (attraction). This property together with the inversion 
symmetry implies that in the stripe \ fa,\ < \U\ — At of the (fa,, fa) -plane it is 
enough to determine the phase diagram at the half-line fa — 0, fa < (or 
fa > 0). 
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Our aim in this paper is to investigate the ground-state phase diagram off the 
symmetry point. According to the state of the art, this is feasible only in the 
strong-coupling regime. From now on, we shall consider exclusively the case 
of the strong positive coupling, i.e. U/t > 4. In this case, it is convenient to 
express all the parameters of H(p e ,pi), in the units of \U\, i.e. we change 
t — > t/\U\, e — > e/\U\, etc, but keep the previous notation. Then, H assumes 
the form 



-* E (■ 

{x,y)i 



C x Cy 



h.c. 



- E 

16 ^ 



(5) 



and the strong-coupling regime is given by the condition t < 1/4. The possi- 
bility of studying the phase diagram in the strong-coupling regime stems from 
the fact that in the open stripe \p e \ < 1 — At, the ion density pi of a ground- 
state configuration determines the electron density: p e = pi for U < 0, and 
p e = 1 — pi for U > 0. Then, one can formally expand the ground-state energy 
Es (p e , Pi) in powers of t, what results in the so called effective Hamiltonian: 



E S (p e , pi) = ~ {pi - p e ) S x ~ 7, + Ate + 1) |A| + 



4 

4 



9^ 
16 



^ s x s y + 
- 1 {x,y)i 



3T _ e_ 
16 ~ 16 



S^Sy + 



? E ^ + ^E(i + 5 S p) + ^ (4) , 

8 (X,„>3 16 P 



(6) 



up to a term independent of the ion configuration and the chemical potentials. 
In (6), P denotes the (2 x 2)-plaquette of the square lattice A, sp stands for 
the product of pseudo-spins assigned to the corners of P, and the remainder 
PS A \ which is independent of the chemical potentials and e, collects those 
terms of the expansion that are proportional to t 2m , with m = 3,4, . . .. The 
expansion (6) (with e — 0) was introduced and the phase diagram, according 
to the expansion truncated at the fourth order, was studied in [10]. It turns 
out however that, in the strong-coupling regime the expansion (6) is absolutely 
convergent, uniformly in A [11,12]. Due to this fact, it is possible to establish 
rigorously a part of the phase diagram (that is the ground states are deter- 
mined everywhere in the (p e , //j)-plane, except some regions), by determining 
the phase diagram of the expansion truncated at the order k, that is according 
to the effective Hamiltonian (p e ,pi). 
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2.2 Construction of the phase diagram up to fourth order. The case of hop- 
ping fermions 



To construct the phase diagram according to the effective Hamiltonian Eg (/i e , /ij), 
we use the m-potential method introduced in [13] , with technical developments 
given in [10,11,12]. By virtue of the comments in the preceding subsection, in 
what follows we set /i e — and fa = fa In the order zero the effective Hamil- 
tonian reads: 



4 0} (o, ^) = -|E( s - + 1 )-^ £ s * s v ( ? ) 

x (x,y)2 

=E^ 0) , 
p 

where 

^ 0) = -f E («* + !)- ^ E ( 8 ) 

Except the point = 0, e — 0, where all the configurations have the same 
energy, the plaquette potentials H^ 1 are minimized by the restrictions to P of 
a few periodic configurations on A. For /i < (/i > 0) it is the ferromagnetic 
configuration 5_ (5V), where s x = —1 at every site (s x = +1 at every site), 
while at the half-line /i = 0, e > 0, besides the ferromagnetic configurations 
S- and 5+, these are the two antiferromagnetic (or chessboard) configurations 
S^ b , where s x = e x , and S° b , where s x = —e x , with e x — 1 if x belongs to the 
even sublattice of A and e x = — 1 otherwise. Moreover, out of the restrictions 
S-\p, S + \p, S^ P , and S° b \ P only four ground-state configurations can be built, 
which coincide with the four configurations named above. Clearly, this is due 
to the n.n.n. ferromagnetic interaction. The phase diagram according to the 
effective Hamiltonian E s °^ (0,/i) is shown in Fig. 1. 

In the next order, which takes into the account interactions proportional to t 2 , 
a n.n. antiferromagnetic interaction appears. The effective Hamiltonian reads 



4 2) (o, u) = -\ E (** + 1) + j E °*°v - YQ E °**v (9) 



=E^ 

p 



(x,y)i (x,y) 2 

(2) 



where, 



HP = -f E (*» + 1) + J E s ^ - 7^ E ***y (10) 

By determining the plaquette configurations that minimize the potentials H p 
one finds the ground-state phase diagram shown in Fig. 2. The boundary 
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Fig. 1. Ground-state phase diagram of Fig. 2. Ground-state phase diagram of 

4 o) (o,/i). 4 2) (°^)- 

between S + and 5_, determined in the order zero, i.e. the half-line /i — 0, 
e > 0, is replaced by the stripe of the width At 2 , centered at this half-line. 
Outside this stripe the ground state remains ferromagnetic. Inside the stripe, 
including its boundary on the line e = 0, there are only two ground-state 
configurations, the antiferromagnetic ones. At the boundary of the stripe, 
given by /i = — 2t 2 ,e > (fj, — 2t 2 ,e > 0), there are exactly three ground- 
state configurations S^ b , and S° b (S + , S^ b , and S° b ). Only at the points (// = 
±2t 2 ,£ = 0) the number of ground-state configurations grows exponentially 
with the size of the underlying lattice (like exp(const|A|)). The Hamiltonian 
i£g 2 ^(0,/i) is known in the physics literature as the Fisher-stabilized Ising 
antiferromagnet [14]. 

To proceed further and analyze the effect of the fourth-order interactions, we 
make use of the inversion symmetry that enables us to fix, without any loss 
of generality, the sign of the chemical potential. From now on, we restrict our 
investigations to the case /x < 0. 

Let us emphasize that we are interested in the phase diagram for sufficiently 
small t. Therefore clearly, it is the neighborhood of radius 0(t 4 ) of the point 
\x = ±2t 2 ,e = 0, in the (//, £~)-plane, where the effect of the fourth order 
interactions can be most significant. In this neighborhood, it is convenient to 
introduce the new coordinates, 5, e, as follows: 

H = -2t 2 + 8t\ £ = et 4 . (11) 

In terms of S, e, the effective Hamiltonian up to the order four can be written 
in the form: 
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E s ] (0, *) = jE ( s * s y + s x + s y + l) 



t 2 



2 



E4 } 



(12) 



t q 



t 2 E l + ^e?, 



with 



■ E + ( 13 ) 



where we have omitted the term t 2 |A|/2. In (12,13), T stands for a (3 x 3)- 
plaquette (later on called the T-plaquette) of the square lattice, whose sites 
are labeled from the left to the right, starting at the bottom left corner, so 
that s 5 is the pseudo-spin of the central site and its left-neighbor pseudo-spin 
is s 4 . The above form of (0, 5) shows manifestly that, in a neighborhood 
of the point 5 = 0, e = whose radius is 0(1), those configurations that 
contain pairs of n.n. sites, (x,y)i, with pseudo-spins s x = s y = 1 cannot 
be the ground state configurations (their energy is larger than the energy of 
the configurations that do not contain such pairs by a large energy of the 
order 0(t 2 )). Therefore, in the considered region of the phase diagram, all the 
effects due to the interactions up to fourth order are described by the effective 
Hamiltonian H^ s \ and the only admissible configurations are those that do 
not contain any pairs of n.n. sites with pseudo-spins taking the value 1. 

To determine the phase diagram due to we use, as before, the m-potential 
method. However, in contrast to the lower-order cases the analysis is not that 
straightforward. The potentials cannot serve as the m-potentials in the 
whole (5, e)-plane. This difficulty can be overcome by introducing the so called 
zero-potentials [10,11,12]. Here, we follow closely [11,12] and introduce the 
zero-potentials such that they are invariant with respect to the symme- 
tries of H and satisfy the condition 

E4'=0. (14) 

T 

Moreover, the zero-potentials are chosen in the following form 

4 4) = E«4 J) , (is) 
i=i 

where CKj are real linear functions of 5 and e (to be determined in the process 
of constructing the phase diagram), while the potentials kj) are given by 



H, 



(4) 



-5 (s 5 + 1) - — E s * s v + 



3 

32 



Y2 E s * s v + ^E( 5s -p + x )> 
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Fig. 3. All the admissible T-plaquette configurations (up to symmetries). 



kfi = Si + s 3 + s 7 + Sg 



k 



,(2) 
(3) 



S2 + S 4 + Sq + Sg 



4s 5 , 
4s 5 , 



: S\S2 + + S3S6 + SgSg + SgSg + SjSg + S4S7 + S1S4 

— 2S2S5 — 2s 5 s 6 — 2s 5 sg — 2s 4 s 5 , 



,(4) 



k T ' — S1S5 + + S5S9 + S5S7 — S2S4 — S4S8 — s$Sq — 

S1S7 — 2s4Sg — 2s2>Sg. 



, (5) _ . . 

K T — S1S3 + S3S9 + S7S9 



(16) 



We note that the potentials kj) are invariant with respect to the symmetries 
of the Hamiltonian H and they satisfy the condition (14). Consequently, the 



same properties are shared by the potentials K T 
as follows 



(4) 



Now, we can rewrite H, 



(4) 
eff 



(4) 



(17) 



T 



and the problem is to find, at each point (5, e), the T-plaquette configurations 
that minimize the potential H^f 1 + K^f 1 . Taking into account the symmetries 
of the Hamiltonian, there are 20 admissible T-plaquette configurations whose 
energies have to be compared. These configurations are displayed and labelled 
(as in [11]) in Fig. 3. The explicit analysis shows that the (5, e)-plane decom- 
poses into five open domains 5_, <S C &, Si, S 2 , and «S 3 , whose boundaries consist 
of straight-line segments, see Fig. 4. We shall adopt the same symbols to de- 
note also the sets of ground-state configurations in these domains. In each 
point (5, e) of a domain, it is possible to choose the values of the coefficients 
a,; such that the set of T-plaquette configurations minimizing the potential 



H. 



(4) 



+ denoted Smt, etc., is the same. Remarkably, out of these min- 



imizing T-plaquette configurations (in a domain) one can construct exactly 
one, up to the symmetries of H , ground-state configuration on the lattice A. 
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S {17, 20} 

CD 




Fig. 4. Ground-state phase diagram of . By the crossing points of boundary-line 
segments we give their coordinates. The numbers in curly brackets, displayed by 
the arrows pointing towards the boundary-line segments, identify the additional 
minimizing T-plaquette configurations. For more comments see the text. 



The obtained configurations turn out to be periodic. Specifically, one finds that 
the minimizing T-plaquette configurations are: S_|t — {!}, <Scb\T — {17,20}, 
S\\t = {2,8}, S 2 \t = {2,5,11}, and S31T = {12,13}, where the numbers in 
the curly brackets stand for the labels, according to Fig. 3, of the minimizing 
T-plaquette configurations, and the T-plaquette configurations that can be 
obtained from the named ones, by applying symmetries of H , are omitted. 
Consequently, the sets of ground-state configurations in the five open domains 
are given by: S_ = {£_}, <S cfe = {S e cb , S° cb }, Si = {S{} 3 ^ , S 2 = {S ] 2 }' J ~X, and 
S3 = {SiYjZi, where the representative configurations of <S C &, «Si, S 2 , and S3, 
are displayed in Fig. 5. The coefficients are given in the Table 1 included 
in Appendix. 



For e = the obtained phase diagram coincides with the phase diagram of 
the spinless Falicov-Kimball model [10,12]. If e < 12, when 5 is increased, it 
crosses some "critical" values, S cr (e), that limit the ^-extent of the domains 
S-, Si, S 2 , S3, and S c b. In particular, S cr (0) = —1/2, 2, 6, 21/2, respectively. 
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S, (ppl/5) S 2 (ppl/4) 




5, (ppl/3) S rt ( Pi =l/2) 

Fig. 5. Ground-state configurations, which appear for /i < 0. 

Similarly, as e is increased, it passes through three "critical" values e cr = 4, 8, 
12, that limit the e-extent of the domains S2, Si, S3, respectively. Thus, there 
is no 1S2 ground-state configurations for e > 4, etc. Above e = 12 the phase 
diagram is independent of e. 

At the boundaries of the domains the situation is more involved. Let S, S' be 
two domains of the considered phase diagram. The set of the minimizing T- 
plaquette configurations at the boundary between S and S' consists always of 
<S|tU<S'|t, but it may contain also some additional T-plaquette configurations 
of minimal energy. If it is the case, a great many of ground-state configurations 
exists at the boundary. Except the boundary between <S_ and S cb , the number 
of ground-state configurations that can be built out of the minimizing T- 
plaquette configurations grows indefinitely with the size of the lattice. 

Specifically, at the boundary between 5_ and S c b the set of the minimizing 
T-plaquette configurations is S_\t U S c b\T- Clearly, only the three periodic 
configurations, 5*_, S* b , and S° b , can be built out of them. 

Then, at the boundary between «S_ and 1S3 the set of the minimizing T- 
plaquette configurations is iS_| T U53| T U{4, 7, 9}. A simple reasoning shows that 
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in any ground-state configuration all the lattice lines with one of the slopes 1 
or —1 are ordered ferromagnetically. On any horizontal line of sites, every two 
consecutive sites with pseudo-spins 1 are separated by at least two sites with 
pseudo-spins —1. That is, the ground-state configurations at the boundary 
"interpolate" between the configurations S3 and <S_. The number of such con- 
figurations grows as exp (const y / ^|A|). The last segment of the boundary of S- 
is the boundary between 5_ and <Si. The set of the minimizing T-plaquette 
configurations is S_\ T U <Si|x U {3,4,9}. Therefore, flipping in Si any set of 
pseudo-spins 1, we obtain a ground-state configuration. Their number grows 
as exp (const I A|). 



The situation at the boundary between <Si and S3 is more intricate. The set of 
the minimizing T-plaquette configurations is Si\tUS 3 \t U {6}. In any ground- 
state configuration the lattice lines with one of the slopes 2, 1/2, —2, or —1/2 
are ordered ferromagnetically. Along such a line, it is possible to build a stra- 
tum, of any width, that consists of elementary cells of Si kind (v^5 x y/E 
squares) (see Fig. 5), then a stratum that consits of elementary cells of S3 kind 
(diamonds), and so on. The number of such periodic configurations grows as 
exp (const yjAj). Moreover, it is possible to build quasi-periodic tilings com- 
posed of square elementary cells and diamonds [15]. 



The three boundaries that remain to be described, have already been discussed 
in the literature. At the boundary between S3 and <S C &, the set of the mini- 
mizing T-plaquette configurations is Ss\t U S c b\r U {19}. In any ground-state 
configuration the lattice lines with slope 1 or —1 are ordered ferromagneti- 
cally. On every horizontal and every vertical line of sites, any two consecutive 
sites with pseudo-spins 1 are separated by either one or two sites with pseudo- 
spins equal to —1 [11]. Thus, the number of such configurations grows as 
exp (const J | A|). 



Next, at the boundary between S2 and £3 the set of the minimizing T-plaquette 
configurations is iS2|tU<S3|tU{6}. In any ground-state configuration the lattice 
lines with one of the slopes 2, 1/2, —2, or —1/2 are ordered ferromagnetically. 
If the slope is ±2, on every vertical line of sites, any two consecutive sites with 
pseudo-spins 1 are separated by either two or three sites with pseudo-spins — 1. 
If the slope is ±1/2, on every horizontal line any two consecutive sites with 
pseudo-spins 1 are separated by either two or three sites with pseudo-spins 
— 1 [11]. The number of such configurations grows as exp(consW|A|). 



Finally at the boundary between Si and S2 , the set of the minimizing T- 
plaquette configurations is Si\t U S 2 \t- Here the situation is analogous to that 
of the boundary between S\ and S3 [15]. 
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2. 3 Construction of the phase diagram up to fourth order. The case of hop- 
ping hard-core bosons 



In the case of hard-core bosons, the creation and annihilation operators at 
site x: c+, c x , respectively, satisfy the anticommutation relations, as fermions, 
but in contrast to fermions they commute at different sites. Consequently, 
with the system of hard-core bosons, described by the Hamiltonian (5), we 
cannot associate a one-particle Hamiltonian and derive the small t effective 
interaction as in the case of fermions. The effective interaction up to fourth 
order has been derived in [12] by means of a closed-loop expansion, and reads: 



E s (fj, e , fr) = ~ (Hi - fl e ) S x - ]- (fit +H e + 1) |A| + 
" x 



^ S xSy + 

{x,y)2 

(x,y)s p 



~t 2 5t A ~ 


E S xSy + 

(x,y)i 


"5t 4 


e 


4 16 


16 


16 



The above expansion is convergent absolutely if t < 1/16 and \fj, e \ < 1 — lQt. 
The latter condition implies that N e + Ni — |A|. Setting in (18) /i e = 0, \Xi = 
—2t 2 + 5t 4 , and dropping the remainder we get the effective interaction 
up to fourth order, (0,5), and the effective Hamiltonian, , defined as 
in (12). Then, repeating all the steps of the construction of the phase diagram 
up to fourth order, carried out for fermions in the previous section, leads to 
the phase diagram displayed in Fig. 6. The diagram in Fig. 6 is described in 
the same way as the diagram in Fig. 4. If the phase diagram due to (0, S) 
is limited to positive e, then it differs significantly from its counterpart for 
fermions, the ^-domain is missing. But the analysis extended to negative 
values of e reveals the iS 2 -domain, and shows that the phase diagram of hard- 
core bosons, limited to e ^ —4, after translating by the vector (2,4) becomes 
identical with the phase diagram of fermions, limited to e ^ 0. 



3 Discussion of the phase diagram and summary 



The ground-state phases of the whole electron-ion system under consideration 
can be distinguished by their ground-state arrangements of ions. Thus, we can 
speak of iS C 6-phase, where the arrangement of the ions is given by the chess- 
board configurations , etc. Let us ignore for the moment the remainder 
in the expansion (6). Then, it follows from the grand-canonical phase diagram 
due to the effective interaction truncated at fourth order that if pi — p e — 1/2, 
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s 



{4, 7, 9} 




Fig. 6. Ground-state phase diagram of H^ s . See description of Fig. 4. 

then the ground-state phase is the <S c &-phase. More generally, consider a se- 
quence of rational ion densities pi = 0, 1/5, 1/4, 1/3, 1/2 (and the correspond- 
ing electron densities p e — 1 — p^. Then, for some densities of this sequence, 
depending on e, the ground-state phases of the system are characterized by 
crystalline, i.e. exibiting positional long-range order, arrangements of the ions, 
given by the ion configurations S-, Si, S2, S3 and S c b, of the corresponding 
density. For e > 12, we can have only two crystalline phases: iS_-phase and 
iScb-phase. At the half- line, given by 5 = 9/2 and e > 12, these two phases 
and only these two ones coexist, i.e. any S that is not in 5_ U S c b has higher 
energy. Since, 
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where i? eff , Hf\ and Ky are evaluated at 5 = 9/2, and E { T ' c j stands for 
the common value of + K^f 1 obtained for S G <S_ U S c b and 5 = 9/2, the 
coexistence of <S_-phase and 5 c b-phase follows from the following inequality: 



T (x,y)i 

C( 4 )t 2 |T_, cb (5)|, (20) 

for any configuration 5 and e > 12. In (20), ifj 4 ^ and are also evaluated 
at 5 = 9/2, T_ c b(S) stands for the set of the T-plaquettes such that S\t 
does not coincide with T-plaquette configurations labelled by 1, 17, 20 (we 
call T- tC b(S) the set of T-plaquette excitations relative to the configurations 
in S- U S c b), \T-,d>(S)\ is the cardinality of T-^S), and is a constant 
(i.e. is independent of S, t, e, and A). The inequality (20), together with 
the variational argument constructed in [16], can be used to prove that for 
any density < pi < 1/2, the ground-state phase is a mixture (a state of 
phase separation) of the 5_-phase and 5 c b-phase. Specifically, the variational 
argument shows that for e > 12 there is a constant such that in any ground- 
state the number |71 jC j,(<S')| of T-plaquette excitations is bounded from above 
as follows: 

c (4) r 2 

|T_, cb (S)K -J|A|. (21) 



£-12 

This means that, for < pi < 1/2 and sufficiently large lattice, any ground- 
state configuration of ions consists of connected regions whose ion configura- 
tions are restrictions to these regions of one of the configurations S-, S^ b , S° b . 
These connected regions are separated by domain walls that involve negligible, 
i.e. 0(|A| -1 / 2 ) (for large lattice A), fraction of sites. 

Now, let us fix the ion density, pi = 1/5 or 1/4 or 1/3, and a sufficiently small 
t. If we choose an initial value of e so that e = e/t 4 is greater than 12 (i.e. 
initially the system is in the state of phase-separation), then by decreasing 
e sufficiently the system is driven into a crystalline state. But obviously, the 
same result can be arrived at, if we fix the value of e and an (sufficiently small) 
initial value of to so that e — eft 4 , — 13, for instance, and then increase t. If not 
the fact that we have ignored the remainder in the effective interaction, 
this would be the announced in the Introduction statement. Thus, it remains 
to be proven that the remainder does not destroy the picture obtained 
with the truncated effective interaction 

First, we shall demonstrate that if we take into account then there is a 
(sufficiently small) constant t such that for t < t the sets of ground-state 
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configurations of ions for densities pi = 1/5, 1/4, 1/3 coincide with the sets 
<Si, S2, and S3, respectively. Specifically, we shall show that for t < t there are 
nonempty two-dimensional open domains S^°, S%°, and S^° that are contained 
in the domains <Si, <S 2 , and S 3 , respectively, and such that in S°°, j = 1,2,3, 
the set of ground-state configurations is Sj. To achieve this we shall construct 
a lower bound for the energy difference E s (0, -2t 2 + 5t 4 )-E Sj (0, -2t 2 + 5t 4 ), 
with Sj G Sj. Clearly, we can restrict the set of all configurations S to those 
that do not contain any pairs (x,y)i with s x = s y = 1. With this restriction, 



E s (0, -2t 2 + 5t 4 ) - E Sj (0, -It 2 + 5t 4 ) = 



I E \{h? + k?)\ s -(hP + kP) 



Ter 3 (s) 



(22) 



Clearly, 



E 

TG75(S) 



H, 



(4) 



T 



if. 



(4) 



(4) 



T 



if, 



(4) 



>Tj\Tj{S)\ 



(23) 



where 7^ (5") is the set of T-plaquette excitations relative to the ground-state 
configurations from Sj and 



Tj = mm 



+ K, 



(4) 
T 



: 5 



(24) 



is a function of (<5, e) on the domain 5^. Following [12] we have also the upper 
bound: 



-i? (4) 



^t 6 Tj\Tj(S)\ 



(25) 



for some constant Tj. From (22,23,25) we obtain 

E s (0, -It 2 + 5t 4 ) - E s . (0, -2t 2 + 5t A ) > t 4 (r, - r-jt 2 ) \Tj(S) \ . (26) 

Note that Sj is an open convex domain in the (5, e)-plane. For (S,e) G Sj, 
the function (<5, e) — > Tj(S,e), being the minimum of a finite set of positive 
linear functions, is concave, piecewise linear and positive. It vanishes only at 
the boundary of Sj. Therefore, there is a (sufficiently small) to such that for 
t < to the value of rjt 2 is less than the maximum of Tj over Sj. Then, the 
level set {(S,e) : Tj(S,e) > rjt 2 } is a nonempty convex open set, contained in 
domain Sj. At this level set the ground-state configurations coincide with the 
configurations from Sj. 

Second, we shall show that for a sufficiently small t the inequality (21) is 
modified by only slightly, so that the conclusions drawn from (21) remain 
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valid. The remainder can be estimated in the following way [16]: 



ij(4) _fl(4) 
5 



S c b 



^ 6 r_, c6 |T_, c& (S)|, 



(27) 



for some constant r_ )C ;,. Therefore, taking into account the remainder i?( 4 \ we 
replace the inequality (20) by 



64 
t 4 



[ £ -(12 + * 2 r_, cb )]|T_, cb (S)|^ 



E (4 4) + 4 4) - 4"' c6) ) + 



t 2 X! i + i? {4) - i? (4) 



S C 6 



^tt 2 |71 c6 (S)|, 



(28) 



for some constant C. As a result, the variational argument of [16] implies that 
there is a constant c such that the following counterpart of (21) holds: 



\T-AS)\ < 



cr 



£ - (12 + 64r_ iCfe t 2 ) 



(29) 



Consequently, for any density < pi < 1/2 and n.n.n. interaction e > 12 + 
64r_ )Cb i 2 the ground-state phase is a mixture of the 5_-phase and iS cb -phase. 

Summing up, let the ion density be one of pi = 1/5, 1/4, and 1/3, and let the 
corresponding electron density be p e = 1 — p^. We restrict our considerations 
to the case pi < 1/2, since the situation for pi > 1/2 can be obtained by means 
of symmetries (see subsection 2.1). Let the hopping intensity t be sufficiently 
small. Then, we have proved that the considered system can be driven from the 
state of phase separation (the mixture of <S_- and iS cb -phases) to a crystalline 
state with a long-range order, which is the iSj-phase whose particle densities 
are the chosen ones, by increasing quantum fluctuations (measured by t) due 
to the hopping electrons. 

Finally, it is easy to see that if the considered model is augmented by a n.n. 
repulsion between the ions, then the above analysis of the phase diagram can 
be carried out without any essential modifications. Thus, we arrive at the 
same conclusions concerning the transition from a state of phase separation 
to a state with a long-range order due to increasing hopping. 

A similar discussion can be carried out in the case of hard-core bosons, with 
one exception. Our analysis shows that the transition from a state of phase 
separation to a state with a long-range order due to an increasing hopping 
occurs for densities pi = 1/5, 1/3. 
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Appendix 

Here we provide the tables of the coefficients ccj, % = 1, ... ,5, of the zero- 
potentials in the case of hopping fermions and hopping hard-core bosons. 

Table 1 



Zero-potential coefficients aj for fermions (e ^ 0) 







«2 


a 3 


a 4 


"5 


S- , e < 8 





(S 
8 


ij-l 

32 T 64 





1 

48 


S_ , e ^ 8 





8 


6 , s 3 
32 " r 128 64 





1 

48 


Si 


34 , e 3 

40 ' 160 80 


4 E , 1 
10 80 + 80 





6 , e 1 
40 + 80 80 


1 

48 


s 2 


74 £ 5 
128 256 64 


94 , E , 3 
64 T 128 " t " 32 





4 £ 5 
128 256 64 


4 £ 7 
128 256 192 


S3, s «: 4 


114 , 5 

144 144 ' 96 


4 
8 





4 | e 7 
144 ^ 288 96 


4 e | 5 
144 288 ' 96 


S3, 4 ^ e 8 


S B 1 1 
12 384 ' 16 


4 
8 





£ 1 
128 16 


" 128 ~^ 24 


S3, e ^ 8 


~i~2 — 96 8 


¥ 








1 

48 


Scb 





4 

8 


1854 , s , 169 
5856 ~r 5856 ^ 3904 





1 

48 



Table 2 



Zero-potential coefficients a« for hard-core bosons (e ^ —4) 





«1 


«2 


a 3 


a 4 


a 5 


S_, e ^4 





4 
8 


S 1 
32 192 





1 

48 


S_,0 4 





4 
8 


4 , £ 7 
32 ~r" 128 192 





1 

48 


Si 


34 , e , 1 
40 ' 160 ' 240 


4 £ i 1 
10 80 ' 80 





4 , £ 1 
40 "•" 80 80 


1 

48 


s 2 


74 £ 7 
128 256 192 


94 | £ | 3 
64 1 128 1 32 





4 5 
128 256 64 


4 7 
128 256 192 


s 3 , e s: 


114 , 11 

144 144 ' 288 


4 
8 





4 | e 13 
144 ^ 288 288 


4 £ | 7 
144 288 288 


S 3 , sCe ^4 


4 e , 5 
12 384 " t " 96 


4 
8 





E 1 
128 32 


— S- + — 

128 ' 96 


S 3 , e > 4 


4 £ i 1 

12 96 12 


4 
8 








1 

48 


S C 6 





4 
8 


1854 , s , 93 
5856 ^ 5856 ^ 3904 





1 

48 
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